Pulsed squeezed-light generation in a waveguide with second-subharmonic generation 

and periodic corrugation 
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Quantum pulsed second-subharmonic generation in a planar waveguide with a small periodic 
corrugation at the surface is studied. Back-scattering of the interacting fields on the corrugation 
enhances the nonlinear interaction giving larger values of squeezing. The problem of back-scattering 
is treated by perturbation theory, using the Fourier transform for non-dispersion propagation, and 
by numerical approach in the general case. Optimum spectral modes for squeezed-light generation 
are found using the Bloch-Messiah reduction. Improvement in squeezing and increase of numbers of 
generated photons are quantified for the corrugation resonating with the fundamental and second- 
subharmonic field. Splitting of the generated pulse by the corrugation is predicted. 
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I. INTRODUCTION 



The process of second-subharmonic generation which 
is just the inverse process to that of second-harmonic gen- 
eration (l[ is interesting not only as a means of frequency 
conversion. It can also serve as an efficient source of 
squeezed light @ which amplitude fluctuations are sup- 
pressed below the limit given by quantum uncertainty 
relations (for a review, see, e.g., 0t9)- Such nonclassi- 
cal light can be emitted both in the fundamental as well 
as second-subharmonic frequency (SSF) fields [7]. Also 
light with nonclassical photon-number statistics can be 
obtained in this process under suitable conditions 0, Q . 
The process of second-subharmonic generation belongs to 
the whole family of optical parametric processes which 
share many common features [hMi~2| |. Among them, 
spontaneous parametric down-conversion with its ability 
to generate entangled photon pairs plays an important 
role 0. 

In homogeneous nonlinear media, the perfectly phase 
matched nonlinear interaction gives the largest possi- 
ble squeezing of amplitude fluctuations. The principal 
squeeze variance of the SSF field asymptotically reaches 
zero for large values of the gain of nonlinear interaction 
[l3l ] . That is why enhancement of the effective nonlinear- 
ity is important as it immediately results in larger squeez- 
ing and lower pumping intensities needed for reaching a 
given value of squeezing. For this purpose, configurations 
with cavities filled by a nonlinear medium have been usu- 
ally used to generate squeezed light (e.g., @,[Hj]). ^ so 
nonlinear waveguides, for which strong spatial localiza- 
tion of optical fields in the transverse plane is charac- 
teristic, profit from enhancement of the effective non- 
linearity. This enhancement has been widely exploited 
when generating photon pairs in classical [i"B4l7j , multi- 
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layer [l8j] , Bragg-reflection [Ti| [2(| and photonic- wire [2l| 
waveguides. Effective nonlinearity in a waveguide can 
also be increased by the use of mode coupling through 
evanescent waves with fields in a neighboring waveguide. 
This occurs due to an addition al sp atial modulation in- 
duced by energy exchange [HI, Hif- Additional spatial 
modulation can also be introduced in a simpler geome- 
try using a linear periodic corrugation at or below the 
waveguide surface [24], HE]. Back-scattering occurring 
at the corrugation modifies electric-field amplitudes of 
the nonlinearly interacting fields which results in the en- 
hancement of effective nonlinearity under suitable condi- 
tions [26 28]. Back-scattering on a periodic corrugation 
has already been exploited for cw squeezed-light gener- 
ation both in the process of second-harmonic [29( and 
second-subharmonic generation pel ] . Scattering on a pe- 
riodic corrugation in a waveguide can also be used to 
enhance second-harmonic generation in Cerenkov config- 
uration [30-32]. Squeezed- light generation in nonlinear 
photonic structures has been discussed in general in [33[ . 

Photonic structures modify in general phase match- 
ing conditions of the nonlinear interaction. Propagation 
constants in waveguiding structures represent a typical 
example. Especially Bragg-reflection waveguides offer 
wide possibilities in this direction [34|, [35|. Also scat- 
tering on a periodic corrugation gives an additional term 
to nonlinear phase- matching conditions [36j-(38|. For this 
reason we need a tool that allows us to reach nonlin- 
ear phase matching conditions for an arbitrary photonic 
structure. Periodic poling 39-42] of x^ 2 ' susceptibility 
has occurred to be extraordinarily useful here and has 
resulted in the so-called quasi-phase-matched nonlinear 
interactions. Using this method even spectrally broad- 
band two-mode nonlinear interaction with femtosecond 
pulses has become possible (43J. 

Due to temporal energy concentration, the pulsed 
regime of the nonlinear process allows to use lower 
pumping powers to observe the needed level of squeez- 
ing. It also brings into attention new features of the 
generated light, namely its spectral modal structure. 
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In the pulsed regime and travelling-wave configuration, 
squeezed-light generation in the considered nonlinear 
interaction has been studied with the help of phase- 
space quasi-distributions or the corresponding Langevin 
stochastic equations |44j. A local oscillator in the 
form of an ultrashort optical pulse is needed to observe 
pulsed squeezing experimentally using homodyne detec- 
tion. The effort to observe the largest possible values 
of squeezing has raised the question about an optimum 
shape of the local-oscillator field [45l - l47j . The Bloch- 
Messiah reduction of the evolution matrices (operators) 
has been proposed for the solution of the modal structure 
(48j . This method has been elaborated in detail in [49ll50j 
for degenerate parametric down-conversion in a BBO 
crystal. Also the relation between the Bloch-Messiah 
reduction and Schmidt decomposition of a two-photon 
spectral amplitude characterizing spontaneous paramet- 
ric down-conversion has been found. The results ob- 
tained in single-pass geometry have been generalized to 
the nonlinear interaction in a cavity [5l[ . Mode structure 
of nonclassical states arising from squeezed states after 
post-selection done with on/off detectors has been ana- 
lyzed in [52( 1 . The discussed effects occurring in temporal 
(or spectral) domains have their counterparts in spatial 
domain, in the transverse plane of the beams. Also here 
eigenmodes typical for the nonlinear interaction can be 
revealed [53j . Circular symmetry of the usual optical 
beams, however, results in different types of eigenmodes 
defined in the transverse plane. Despite this, a close sim- 
ilarity in the behavior of fields in the spectral and spatial 
domains can be found. As an important example, cor- 
relations between intensities of the interacting fields in 
their transverse planes can be mentioned [54], [55[ . 

Here, continuing the investigation in [28j | we pay at- 
tention to the pulsed SSF generation in a waveguide 
with linear periodic corrugation on its surface that causes 
back-scattering of the interacting fields. Assuming strong 
pulsed fundamental field, we study squeezed-light gen- 
eration in the SSF field. We pay our attention to the 
enhancement of nonlincarity and the related increase of 
squeezing applying three different approaches. We uti- 
lize a sophisticated perturbation approach, the method 
of Fourier transformation in non-dispersion field propaga- 
tion, and numerical solution for the general case. Squeez- 
ing is characterized by a principal squeeze variance intro- 
duced in 0, H, @ that is determined for suitable spectral 
modes. 

The paper is organized as follows. In Sec. II, a multi- 
mode quantum model of the nonlinear interaction us- 
ing the appropriate momentum operator and the related 
Heisenberg equations is presented. Sec. Ill is devoted 
to the perturbation solution of the model that is divided 
into three parts. A general perturbation solution is found 
in Subsec. Ill A, Gaussian spectral approximation to the 
solution is applied in Subsec. IIIB, and squeezing is an- 
alyzed in Subsec. IIIC. If inter-mode dispersion is omit- 
ted, the model can be solved by the Fourier transform, 
as shown in Sec. IV. Discrete formulation of the model 



is elaborated in Sec. V devoted to the numerical solu- 
tion. Also quantities useful in the characterization of the 
generated field are introduced in this section. Discussion 
of squeezing and appropriate spectral eigenmodes is con- 
tained in Sec. VI. Whereas the model with non-dispersion 
field propagation is analyzed in Subsec. VIA, the results 
of the general approach are studied in Subsec. VIB. Con- 
clusions are drawn in Sec. VII. An optimum mode profile 
giving the maximum pulsed squeezing is found in Ap- 
pendix A. 

II. QUANTUM MULTI-MODE MODEL OF 
SECOND-SUBHARMONIC GENERATION 

We consider a nonlinear waveguide made of LiNb03 
(see Fig. [I} in the configuration that allows to generate 
a second-subharmonic (SSF) field using \^ nonlinearity 
and pumping, e.g., by the second-harmonic of a Nd:YAG 
laser at the wavelength of X s = 1.064 x 10~ 6 m. Under 
a suitable choice of waveguide parameters, the waveg- 
uide is single-mode for both the fundamental and SSF 
field and allows for efficient nonlinear interaction. A lin- 
ear corrugation fabricated at the top of the waveguide 
leads to back-scattering of the interacting fields which 
results in the enhancement of electric- field amplitudes in- 
side the waveguide under suitable conditions. This leads 
to effective increase of the nonlinear interaction and gives 
larger amount of squeezing of the SSF light. Quasi-phase 
matching of the nonlinear interaction is guaranteed by 
periodic poling with an appropriate poling period. A de- 
tailed description of the waveguide was given in (28|. It 
has been shown that the investigated waveguide can be 
described by the following parameters: propagation and 
coupling constants of the fundamental and SSF fields and 
constants characterizing the nonlinear interaction occur- 
ring among both the forward- and backward-propagating 
fields. When losses inside the waveguide caused both by 
absorption and scattering of the light outside the guided 
modes are neglected we can describe the nonlinear inter- 
action in the waveguide as follows. 

According to quantum theory, the electric-field vec- 
tor operator amplitudes E a (x, y, z, t) (a = p, s) at time t 
and spatial point (x, y, z) inside the waveguide can be de- 
composed into harmonic plane waves with mode operator 
amplitudes a a in the Heisenberg picture @, Q : 

/>oo 

E a (x,y,z,t) = / duj a ~E a (x,y, z,u a ) exp(-iuj a t), 
Jo 

(1) 

E a (x, y, z,uj a ) = i[a aF (z,uj a )e a (x,y,LL> a ) 

+ a aB (z,u}a)e a (x,y,oj a ) - H.c] , a = p,s. (2) 

In Eq. ©, a aF (z,uj a ) [a aB ( z i w n)] denotes an annihila- 
tion operator of the mode with frequency u a in field a 
propagating forward [backward] . We note that the mode 
vector functions e a as well as the corresponding propaga- 
tion constants j3 a depend on the frequency uj a and their 
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FIG. 1: Scheme of a periodically-poled nonlinear waveg- 
uide with a rectangular transverse profile (thickness t and 
width Ay) and length L having \^ susceptibility. A lin- 
ear corrugation with period A; and depth ti occurs on the 
waveguide upper surface; A„; denotes period of nonlinear 
poling. The waveguide is made of LiNbOa which optical 
axis coincides with the x axis. Four fields interact inside 
the waveguide: forward-propagating fundamental (electric- 
field amplitude A PF ), backward- propagating fundamental 
(A PB ), forward-propagating second-subharmonic {A SF ), and 
backward-propagating second-subharmonic (A SB ) fields. 



form can be found in [281 ] . 

Evolution of the nonlinearly-interacting quantum op- 
tical fields inside the waveguide is described by the fol- 
lowing momentum operator G 0, [23[ : 



G(z) 



poo 

^2 ^2 duJ a (±l) b hl3 ab (uj a )ai b (z,uj a )a ab (z,uj a ) 

i=p,s b=F,B 

F [ ^ / dcj a hK a (uj a )exp (i^-zj 

a=p,s J V A a J 

x al p (z, u> a )a aB (z, ui a ) + H.c. 

^ 2i duj s / du> s (±l) b K n i A (Lo s ,u s ) 
b=F,B Jo Jo 

x exp {i-j^^j ai b {z,u s )at b {z,uj' s )a Pb (z,uj s + lo' s ) 



H.c. 



(3) 



(±1)_f = 1 and (±1)b = — 1- The linear coupling con- 
stants K p {ujp) [K s (oj s )\ describe scattering on the corru- 
gation with period A p [A s ] that leads to the coupling be- 
tween the fundamental [SSF] fields propagating forward 
and backward. Under suitable conditions, this scattering 
leads to increased values of electric-field amplitudes in- 
side the waveguide. However, this increase is much lower 
than that found in layered photonic band-gap structures 
[56M59I ] . On the other hand, the nonlinear coupling con- 
stants K n i tq (uj s ,uj' s ) characterize the nonlinear interac- 
tion among co-propagating fields and obey the symme- 
try relation K n i tq (u) s ,u s ) = K nl:q (J s ,uj s ). Constant A„ ; 
gives the poling period whereas the integer number q de- 
termines a harmonic frequency employed in the quasi- 
phase-matching. Planck constant is denoted as h and 
symbol H.c. substitutes the Hermitian conjugated term. 
The Heisenberg equations, dX /dz = —i/h[G, X] for an 



arbitrary operator X, written for the momentum opera- 
tor G in Eq. ([3]) attain the form: 



da SF {z,Lo s ) 
dz 



ij3 s (u) s )a SF (z,oj a ) 
+ iK s (uj s )exp (i-^-^j a SB (z,L0 s ) 
+ 4/ duj s K n ^ q (uj s , ui s ) exp ( i— — z 



A 



da SB {z,w s ) 
dz 



da PF (z,uj p ) 



dz 



nl 

x a PF (z,u s +u' s )al F (z,u}' s ), 
-i/3 s (u s )a SB (z,u s ) 

- iK*(uj s )cxp (^i-j^- 2 ^ a SF (z,uj s ) 

f°° 1 1 ( 2irq 

-4/ dto s K n iq(ui Si lo s ) exp [i—, — z 
Jo V A m 

x ci Pb (z,uj s +uj s )dl B (z,uJ s ), 
i(3 p (u p )a pF (z,Lj p ) 



2tt 



+ iK p {uj p ) exp ( i—z ) a PB (z, w p ) 



-2 J duj s K* lq (uj s ,uj p ~uj s )exp^-i^-z^ 
x d SF (z, uj s )a SF (z, u p -oj s ), 

da PB (z,uj p ) 



dz 



-i(3p(u> p )a PB (z,uj p ) 



iK*(u} p )exp [ ~ i A ~ z ) a PF {z,u; p ) 



+ 2/ dujsK^ (uj s ,uj p - uj s )exp [ -i^^-z 

JO V A rd 

x a SB (z,u) s )a SB (z,L) p - u) s ). (4) 

Creation and annihilation operators of the incident 
fields are assumed to fulfil the boson commutation re- 
lations, i.e. 

[a aF (0,oj a ),al, (0,uv)] = S a ^S(uj a - u a ,), 
[d aB (L,u a ),dl^(L,u a ,)] = 5 at a>S(u) a - u) a >), 

a=p,s. (5) 

The remaining commutators are zero. It has been shown 
in |60| | for quadratic momentum operators G that also the 
output operators obey the boson commutation relations: 

[a aF (L,U a ),dl, (L,UJ a >)] = S a .a'S(uJa - U) a >), 

[a aB (0,w a ),a;^(0,a; o ')] = S a>a 'd(u a - u a >), 

a=p,s (6) 

and commutators not mentioned in Eq. ([6]) are zero. 

We note that the nonlinear operator equations written 
in Eq. Q have one integral of motion arising from the 
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conservation of energy flux: 

d_ 

dz 



duj s al F (z, uJ s )d SF (z, lu s ) 
du s al B (z,u)s)a ag (z,u) a ) 

OO 

+ 2 I dujpdl F (z,uj p )d PF (z,ujp) 
dw p d pB (z, u> p )d PB (z, uj p ) 



0. (7) 



It is convenient to introduce new operators A that take 
into account the harmonic spatial evolution induced by 
the corrugation present in general in both the fundamen- 
tal and SSF fields, d aF (z,uj a ) — A aF (z, ui a ) exp(i7T2:/A a ) 
and a aB (z,u: a ) = A aB (z, uj a ) exp(-i7rz/A a ) for a = p,s. 
Equations ((4]) written for the operators A take the form: 

dA SF (z,uj s ) .5 s (u s ) * * 

= i — z — A SF (z, uj s ) + iK s (u; s )A Sb (z,uj s ) 

dz 2 

+ 4 / duj s K nLq (uj s ,uj s )exp[iS n i tq z} 
Jo 

x A Pf (z,uj s +ui 3 )Al p (z,ui s ), 
dA SB (z,ui s ) .S s (uj s ) 



dz 



A SB {z,u s ) - iK*(uj s )A SF {z,uj s ) 



du) s K n i t q(oj s ,u s ) eKp[-id n i tq z] 



(8) 



x Ap B (z, u s + u s )Al B (z, cj s ), 

dA PF (z,ujp) . Spjujp) 2 \,-irr\Xi \ 
j- = i — ^ — A PF (z,ujp) +iKp(u}p)Ap B (z,u}p) 

P OO 

-2 / dui s K* d 9 (w s ,w p - uj s )exp[~i5 nLq z] 
Jo 

x A SF (z, uj s )A Sf (z, u) p -u) s ), 
dA PB (z,ujp) _ .Sp(ujp) 



dz 



A Pb (z,uj p ) - iK*(u}p)A PF (z,ujp) 



+ 2 duisK^ q (uj s ,uip ~ uj s )exp[iS n i iq z] 
Jo 

x A SB (z,u s )A SB (z,L)p - LJ S ). (9) 

Linear phase mismatches 8 p , 6 S and nonlinear phase mis- 
match 8 n i, q (q — 0, ±1) are defined as: 

2- 



<5 a (w a ) = 2j3 a {u) a )-— , a = p,s, 

A-a 

6 — - 2— + — 

Ap A s A n l 



(10) 



We note that whereas q was used for the nonlinear inter- 
action among the forward-propagating fields, — q was cho- 
sen for the interaction among the backward-propagating 
fields in Eqs. ([8]) and ©. We also note that equations 
similar to © and (O can be derived from the wave equa- 
tion considering classical fields and invoking paraxial ap- 
proximation [12j . 



We assume that the fundamental field is strong and its 
depletion due to the interaction with the SSF field can 
be neglected. In this case, the equations in © become 
linear and break into groups with two linear equations 
for a given frequency lu p : 



dA PF (z, uj p ) 
dz 

dA pB (z,U!p) 

dz 



S (lu ) 

i P - P A PF (z,u>p) + iK p (ujp)A 



2 

.Sp(uip) 



A VB (z,u p ) - iK*(uj p )A PF (z,ujp). 

(11) 



The solution of Eqs. can be easily obtained in the 
matrix form 



A PF {z,ujp) 
A PB (z,ujp) 



± 



B p (uj p ) exp[±iA p (w p ) 



A PF (0,ujp) 

Kb (L,U P ) 



(12) 



using the eigenfrequencies ±A p , A p (u> p ) = 
y^pVp)/ 4 ~ I^p( w p)I 2 - Symbol J2± occurring in 
Eq. (1121) means the summation over the terms differing 
in their signs. We note that the boundary condition for 
the backward-propagating field is chosen at z = L. The 
matrices B^ujp) are defined as 



B±{u v ) = Dp 



(± 5 -f + Ap) exp(=RApL) ±K P 



TK;exp(TiA p L) Tf + A p 

(13) 

and D p (ujp) — [2A p cos(ApL) — iS p sin(Api)] 1 . 

If the corrugation is missing, the matrices B p in 
Eq. (|13[) take a simple form obtainable in the limit 
K p — > together with A p — > oo: 



1 






exp(i/3pL) 



(14) 



It also holds that 6 P = 2/3 p , A p — f3 p and D p 
exp(i/3pL)/(2/3p) in this limit. 



III. PERTURBATION SOLUTION 

Perturbation solution of two equations ([5]) can be 
found after substituting the solution for the fundamental 
field contained in Eq. (|12[) . In the perturbation approach, 
the inquired solution A Sb (z,uis) for b = F,B is expressed 

as Xl^Lo {z, uj s ) where an rt-th term A s ™\z, u) s ) is 
proportional to . 
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A. General solution up to the first order in 
nonlinearity 

The zeroth-order terms Ay) (z, uj s ) and A SB (z,u} 8 ) are 
given as a solution to the equations 



It can be again written in the matrix form using the 
matrices Bf(uj s ) from Eq. (JTTJ) : 



A Sf ,(z,u s ) 
A Sb (z,uj s ) 



= ^B±(w a )exp[±iA a (u) a )i 



dA ^ Us) = -i^iW(, )Ws )_^ K )i(0)( 2)a;s) . 

(15) 



" i SF (0,w s ) " 


+ 


F sf (z,oj s ) 


_ A sb (o,uj s ) _ 




. F SB (z,u) s ) _ 



(21) 



This solution can be written in the compact matrix form: 
= Bf(uj s )exp[±iA s (uj s )z} 

(16) 



A SF (z, uj s ) 
. A Sb (z,uj s ) 



A Sf {0,uj s ) 
A sb (0,uj s ) 



A s (w s ) = ^5 2 s (lu s )/A-\K s (lo s )\ 2 . The matrices B±(w a ) 
are expressed as 



1 

2A S 



±4f + A S 



i 2 



(17) 



The equations for the first-order terms As) (z, u) s ) and 
As B (z,u) s ) have a more complex structure: 



dA. 



(i) 



(z^ = ,%)i(i) (z , Ws) + ^ sK) i(i)(,, Ws ) 



dz 



duj s K n i iq (u> s , uj s ) exp[iJ„ iig z] 
x A Pf {z,uj s +J s )A^{z,uj' s ), 



dA s ^(z,u) s ) . S s {uj s ) m) m) 
J z = ~t 2 A sJ( z ^s) - iK s (uj s )A { s J(z,uj s ) 



du} s K n i jq (ui s , u> s ) exp[-i<5„; !(Z z] 

x A PB (z^ s + J s )A^{z,J s ). 



(18) 



Equations (|18p for a fixed frequency cj s represent a cou- 
pled set of two linear differential equations with nonzero 
right-hand sides. They can be formally written as follows 



dA{z) 
dz 



iKA{z) + F(z), 



(19) 



where A and F are vectors with 2 elements and K is a 2 x 
2 matrix. Direct inspection confirms a non-homogeneous 
solution to Eq. (fH?)) in the form 



dz'A hom (z-z')^V) 



(20) 



provided that (L4 hom (z)/(dz) = iKA hom (z). 

The sum of the zeroth- and first-order solutions gives 
us the solution to the equations ((5J) valid up to K n ^ q . 



Operators F Sb (z,oo s ) for b — F, B describe the solution 
proportional to the first power of nonlinear constants 
K n i^ q . They are expressed for z — L below in Eq. (l23l) . 

Now we write the solution in Eq. (l2~Tj) for z = L and 
partially invert the obtained linear relations in order to 
find the input-output relations of the waveguide. The 
result can be expressed as 



A SF (L, lo s ) 

A sb (o,cjs) 



2A S 2iK s sin(A s L 
2iK* sin(A s Z) 2A S 



A SF (0,U) s ) 

A SB (L,u s ) 

F SF (L, u>s) - 2iK s D s sm(A s L)F SB (L, oj s ) ' 
-2A s D s F Sb {L,cj s ) 



(22) 



In Eq. ([2"2]h the definition D s (oj s ) = [2A S cos(A s L) - 
iS s sin(A s .L)] _1 has been used. The expressions on the 
right-hand side of Eq. (|2"2l can be substantially sim- 
plified if the signal field fulfils the resonance condition 
A s w mir/ m = 1,2, . . ., that gives the maximum en- 
hancement of its electric-field amplitudes in the waveg- 
uide. In this case, sin(A s L) m 0. 

The operators F Sb (L,uj s ) at z — L are obtained in the 
form: 



F Sa (L,ij s ) = 
(±1). 



x exp 



(±1)„4 / dw' s K n i iq (w s ,u's 
Jo 

iS n l,qL 



E 

± s b=F,B 



x E E b pm^p 

± P c=F,B 



E ^V x p 

=F,B 

iA p (uj s + lu' s )L 



± 



A pc (u> s + uj' s ) 



x 



E ' ' 

d=F,B 



Bh*, exp 



iA s {uj' s )L 



AlM)s>^{^'s), 

F,B, (23) 



where 

S*"" ± *(w«,wi) = 2Lsinc[[(±l) <J nZ ,,± p A p (w s + fa£) 
± s A s (a; s )± s A s K)]i/2], a = .F,fl, (24) 

sinc(a;) = sin(a;)/x. The matrices B^(u> s ) introduced 
in Eq. (|23|) are defined analogously to those given in 
Eq. (fT3")) and characterizing the fundamental field. We 
note that the expression in Eq. (|23|) contains only one 
half of all possible terms. The missing terms are far from 
the quasi-phase-matching conditions and thus give neg- 
ligible contributions. 
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B. Gaussian spectral approximation in a resonant 
term 

The right-hand-side of Eq. (|23|) giving F Sa (L,ui s ) is 
composed of four terms that differ in the signs of eigen- 
values A p and A s (resolved by the symbols ± p , ± s ). 
Quasi-phase-matching conditions in the waveguide are 
such that they emphasize only one out of the four terms. 
The remaining terms give small contributions. That is 
why we pay attention only to one of them. We also con- 
sider only the incident forward-propagating fundamen- 
tal field in a multi-mode coherent state with a Gaussian 
spectral shape and spectral phase variations such that 
the fundamental field in the middle of the waveguide has 
the same phase along the spectrum: 



A Pf (0,uj p ) = £ p 



2tt 



■ exp 



0\2 



x exp 



iA p (u) p )L 



A pb {L,lo p ) = 0. (25) 

The fundamental pulse has amplitude £ p , duration t p and 
carrying frequency uj~,. 

We further assume that both the fundamental and SSF 
spectra are not too wide and so the propagation constants 
Pa(ua) can be approximated by their second-order Taylor 
expansions: 

AK) = f3 a +/3 la (u a -Lu° a )+/3 2a (u a -Lj a ) 2 , 

0°a = /3a(w2), 



R - l d ^ a 
l la ~ i\ duji 



i = 1,2; a=p,s. (26) 



In this approximation, the eigenvalues A a (w a ) can be 
expressed as 

A a (uj a ) = A° a + A la (uj a - uj° a ) + A 2a (u; a - u° a ) 2 , 
K = A ( W °), 

/3 la (/3° - 7r/A ) 



Ala 

A 2a 



2A° 



2(A° a ) 3 



a = p,s. 



(27) 



Efficient nonlinear interaction occurs if quasi-phase- 
matching conditions for the central frequencies w p and 



,o - ,.,o 



ojp/2 are fulfilled, i.e. 

S n i, q ± p A p {uj p)T s 2A s (w°) =0. 
,± D ,± 



(28) 



The function S a p ' s in Eq. (|2~41 can then be rewritten as 



S^'^^X) = 2£sinc^ 



(zbpAip ± s Ai s ) 



A 2s (uj s 



l\2 



(29) 



assuming Ai a 3> A2 a Aw for a — p, s; Aco characterizes 
the fields' spectral width. 

The use of Gaussian approximation to the sine function 
in Eq. (|29j) [sinc(aa; + fiy 2 ) « exp(-a 2 a; 2 /5 - |/% 2 /3)] 
for constants a, (3] allows to derive the useful relation: 



S ± ^(co s ,uj' s )Ap F (0,uj s +u J / s ) = 2L£ P1 



'2it 



x exp 



iA p (uj s + uj' s )L 



*(«.,«.)■ (30) 



In Eq. (|gQ]l. the function 

(uj s -uj' s ) 2 (uj s + uj' s - w°) 2 



*(w s ,w s ) = exp 



A+ 



1 

aT 
i 

AT 



A 2 ,£ 
12 ' 

(zbpAip ± s Ai s ) 2 i 2 t 2 



20 



4 ' 



(31) 



(32) 



determines the structure of spectral modes. 

These spectral modes can be revealed using the 
Schmidt decomposition [49j,l61] of function $: 



$(w s ,w s ) = ^2 n n <j) n (uj s )<j) n (uj s ). (33) 

n=0 

Eigenvalues fi n introduced in Eq. (j3"3")l take the form: 

2 * 



i/2 



A+A- 



(34) 



A_ 



A. 



Eigenmode spectral functions <f> n can be expressed in 
terms of the Hermite polynomials H n : 



<^n(w s ) 



on f r exp [~ T s(^ 
x U n (T s [uj s -uj]), 



.0\21 



l-6> 2 



(35) 
(36) 



For the considered waveguide, A + <C A_ and thus a typ- 
ical time constant r s of the SSF field defined in Eq. (|3"6")l 
can be approximated as: 



(37) 



The Schmidt decomposition in Eq. (|33|) allows us to 
transform Eqs. (|22[) written for the 'continuous index' 
lu s to those related to spectral modes. The appropriate 
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unitary transformation of the input and output operators 
takes the form: 



duj s (j) n (uj s )a SF ,(0, lu s ) exp 
iA.(u.)L 



x exp 



dw s 4> n (u) s )a SB (L,u s ) exp 
iA s (oj s )L 



~2A~ 



inL 



x exp 



dw s (l) n (LJ s )a sp (L,u> s ) exp 
iA s (w s )A 



2A.< 



«7rA 
~2A7 



x exp 



da; s 0„(a; s )a SB (0,w s ) exp 
iA s (u s )L 



inL 
2A7 



x exp 



(38) 



When weak spectral dependencies of multiplicative fac- 
tors occurring in Eqs. ((22]) and (l23l) are neglected, the 
transformed Eqs. (|22[) valid either close to the resonance 
condition or for non-scattered SSF field are expressed as: 



s F: n 
~Ollt 



2A°£>°exp 



p 

P 

1 SB,n 







' 1 


" 




■ -in 


As 







1 




. a s B ,n _ 



(39) 



Defining suitable nonlinear coupling constants £ n for 
individual spectral modes, 



£„ = 8LK nl , q (u>° s ,u° s ) — ^£ pAt „, 
\ V2w 



(40) 



the operator coefficients A SFi „ and F SB ^ n introduced in 
Eq. ([23) take the form: 



A, 



D°D°* 

s 

2A° 



(5° 



<52 



x ± p f + A£ ± pJ RT-(^r 



-inf 



xif°* ±* +A° ipi^* ^t,A 



5? 



-±- A° 
2 



-i„t 



SB ,» [ ' 



As s ,n £n exp 



i5nl,qL 



exp(± s zA°A)A^|D°| 2 



x (| T S A° 



sint 



f -A°)(i^ 2 



T P — ,, 



52 



Tpif p °* f A° 



SR,I1 f ' 



(41) 



0°=A>(w°). 

If the corrugation is only in the fundamental field, the 
expressions for operator coefficients A SF: „ and A SB ,„ in 
Eq. (jH]) simplify: 



A, 



Cnexp 



D° p exp(-if3° s L) 



6° 
- p ~2 



x ± D ^ + A° \af} 



A, 



-Cn exp 



p I SF,m 
iS n l,qL 



D°exp(-i/3°L) 



x (T ? ^)^, 



p J"SB,n- 



(42) 



On the other hand, the corrugation resonating with 
the SSF field leads to the following expressions 



£ n exp 



iSnLqL 



D 



" 



exp(i#£)^(±.f+A 



2 



_|_ _£. i A I xinf i fr-0*-int 



A, 



-£n exp 



exp(i$°L) exp(± s zA°A) 



«5 



x l^l 2 ^* N Tfl A° ) ai n F t n - *T°*<L } . (43) 



Multiplicative factors occurring in the expressions in 
Eqs. (|4"D — [4"3"|) describe the enhancement of nonlinear in- 
teraction due to scattering of the fundamental and SSF 
fields. This enhancement can be quantified by the expres- 
sion (± a <5°/2 + A°)/(2A°) for a = p,s. Let us consider 
first the corrugation present in the fundamental field. In 
this case, the quasi-phase-matching condition written in 
Eq. ([25]) takes the form 



<52 



rnat,0 p 
°nl,q 



2 ± p a; = o, 



(44) 



where <5°^ = /3° — 2/3" + 2irq/k n i describes the natural 
quasi-phase mismatch. Efficient quasi-phase-matching 
can be reached only if and 6® have the same 

sign [28[. Considering a positive value of phase mis- 
match , the sign + [-] in Eq. (|4"4"]l is suitable for 



€7 < K\ K? q > The enhancement of 

fundamental-field amplitudes is then described by the ex- 
pression 



g/2 + Aj i±„i i > ! 



2A° 



2A° 



(45) 



On the other hand, a negative value of phase mismatch 
SnTq requires the opposite choice of signs in Eq. (|44]) and 
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the enhancement of fundamental-field amplitudes can be 
quantified by the expression 



S° p /2- 



A 

p 



2A° 



rnat,0 
°nl,q 

~2A° 



> 1. 



(46) 



The presence of corrugation in the SSF field needs the 
following quasi-phase-matching conditions: 



(47) 



In this case, the quantities 6™^' and 6® have to differ 
in their signs. The enhancement factors of the SSF-field 
amplitudes can be analyzed similarly as for the funda- 
mental field. 



C. Principal squeeze variance 



The formulas in Eq. (1501) considerably simplify for the 
corrugation in the fundamental field only: 



\ - i _ Km 

A s F ,n 




6° N 



A° 

p 



(51) 



On the other hand, the presence of corrugation only in 
the SSF field leaves us with the expressions 



2(A°) 2 

62 



s2 



x ! ± s f + A° s ±K° s * 



2(A°) 2 



(5° 

±sy + A s ±« ^ s 



(52) 



The enhancement of electric-field amplitudes due to 
the presence of the corrugation results in larger squeezing 
of fluctuations of these amplitudes. The suppression of 
amplitude fluctuations can be quantified by a principal 
squeeze variance A 0, H, [23| that can be determined along 
the relations 



A, 



l + 2(B Sb -\C Sb \) 



B Sb = (Afit Ac 



C 



s b - ((A<) 



(48) 

F, B. (49) 
means the 



In Eqs. (I49I) . Aa = a — (a) and symbol I 
quantum mechanical mean value. 

Substitution of the expressions in Eq. (|4Tj) into 
Eqs. (|48[) and (|4"9")l provides the formulas for principal 
squeeze variances X SF , n and A S£j! „ of individual spectral 
modes, 




62 



(50) 



The enhancement factors discussed for the fundamental 
field in Eqs. (j45|) and (j46|) are clearly recognized in the 
expressions ([5D|1. Also the factors 5°/ 2 ± s A° + i\° and 
S° s /2 ± s A° s + K°* found in Eqs. © significantly con- 
tribute to the enhancement of squeezing of SSF electric- 
field amplitudes due to the same signs in front of 8® and 

K2. 



IV. NON-PERTURBATION SOLUTION FOR 
NON-DISPERSION PROPAGATION 

The linear coupling constants K p , K s and nonlinear 
coupling constants K n i >q depend usually only weakly on 
frequencies in a relatively wide interval. Weak frequency 
dependence of the coupling constants K p , K s and K n i i 
in the considered range [approx. 80 nm (40 nm) for 
the fundamental (SSF) field] is shown in Figs. [5^, 
and [4^,, respectively, for the analyzed waveguide. Also 
inter-mode dispersion both in the fundamental and SSF 
fields can be in the first approximation neglected (see 
Figs. [2b, [3b, and Hp for frequency dependence of the 
linear phase mismatches S p and S s and nonlinear phase 
mismatch S n n). We note that the needed spectral range 
of the SSF field roughly depends inversely proportion- 
ally on the waveguide length L. The solution of opera- 
tor equations considerably simplifies when the frequency 
dependencies are omitted and even certain analytical re- 
sults can be obtained. 

In this case, Eqs. ((8]) and ([9]) for the operator ampli- 
tudes A ab (z,w), a = p, s, b = F,B can be decoupled in 
their frequency 'index' using the Fourier transform: 



duj a A ab (z,u> a )exp(-iu! a T a ), (53) 



Aa b (z,T a ) 



A ai (z,Ua) = 7T I dT aA ab (z,T a )exp(iuj a T a ).(54) 
After the transformation, Eqs. © attain the form: 

J- = i^A Sf {z : t)+iK s A Sb {z 1 t) 

+ ^ K ni, q exp[iS n i tq z]A pp (z, t)AI f (z, t), 

= -i S jA SB (z,T)-iK s *A Sp (z,T) 
~ AK li, q exp[-^„;, g z]i PB (z, t)A\ b (z, t); (55) 
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FIG. 2: (a) Absolute value of linear coupling constant K v and 
(b) linear phase mismatch 5 P — Sp for the fundamental field 
as they depend on relative frequency ui p /u]p\ t — 5 x 1CP 7 m, 
ti = 5 x 10~ 8 m. 



FIG. 3: (a) Absolute value of linear coupling constant K a 
and (b) linear phase mismatch S s — 5® for the SSF field as 
functions of relative frequency uj s /u] p \ values of parameters 
are the same as in Fig. [2] 



K° lq = K nLq (oj^tu°). The solution of Eqs. © writtenin 
Eq. (|12[) is transformed in the considered approximation 
of strong fundamental field as follows: 



Kb 



= ]>> p ±0 exp[± ? A£z] 



A Pf (0,t) 



and 



p 



(w°) defined in Eq. (Jl3 



(56) 

The 
and 



fundamental-field operator amplitudes A PF (0, r) 
A PB (L, r) describe the incident pulses. 

We further discuss the solution to Eqs. ([53)1 for the cor- 
rugation present either in the fundamental or SSF field. 



A. Corrugation in the pump field only 



In this case, two operator equations (1551) are indepen- 
dent. Moreover, we are interested in a solution close to 
the resonance where the nonlinear terms in Eqs. (|55[) give 
considerable contribution. We further pay attention to 
the SSF forward-propagating field and assume only the 
incident fundamental field at z = 0. Using Eq. (|56"|) the 
first equation in (1551) can be rewritten for two different 
resonant conditions indicated by upper indices ±, 



dA SF (z, t) 
dz 



i|i SF (z,r)+4^ ± (r) 



(57) 



The effective nonlinear coupling constant i^ ± (r) = 
K-ni q^pFF^PF (0) T ) incorporates the enhancement of 



b) 



1.04 



1.02 



-5 1.00 



3 



0.96 



1.04 



1.02 



0.98 



0.96 



IKJ^m- 1 ) 




0.96 0.98 1.00 1.02 1.04 

co /oo 

s s 

5, -5 M °(10 6 m" 1 ) 

n .1 n .1 * ' 




0.96 0.98 1.00 1.02 1.04 

oo /to ° 



FIG. 4: Contour plots of (a) absolute value |-Kni,i| of non- 
linear coupling constant and (b) nonlinear phase mismatch 
5 n i,i — $nl l as they depend on relative frequencies uis/^p and 
u} s /uip. Values of parameters are the same as in Fig. [2] 
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nonlinear interaction due to the fundamental-field scat- 
tering. 

The substitution A Sf (z,t) — A sp (z,r) exp[i(<5°; ± 
A° p )z/2], At F (z,r) =^ SF (z,r)eM^(S^ q ±A p )z/2] in 
Eq. (1571) leads to differential equations with constant co- 
efficients obeyed by the operator amplitudes A SF and 
A.\ . Their solution transformed to the original opera- 
tors can be written as: 



A Sf (L,t) = U ff (t)A Sf (0,t) + Vff(t)AI f (0,t), 



(58) 



Uff{t) = i exp(i/3°L) exp (i^j^ 



1 - 



2A±(r) 
2^) 



exp[A ± (r) J L] 
escp[-X ± {T)L] 



V ff (t) = 



exp(ipgL) exp « 



A±(r) V 2 

x [expp^ (t).L] -exp[-A ± (r)L]] 



(59) 



The phase mismatches and eigenvalues A 1 * 1 are given 
along the expressions: 



5°± A 



A ± (r) = Vl6|^ ± (^)l 2 - ^ ±2 /4. 



(60) 



Ideal phase matching leads to fi ± = 0. 

The inverse Fourier transform (l54l) of the expression 
for operator amplitude A SF (L, r) in Eq. (155)) provides 
the spectral operator amplitude A Sf (L,u> s ) that equals 
the operator amplitude a° F (cj s ) in the limit A s — > oo, 



<1uj s Uff(u s ~ uj s y s n F (u s ) 
du' a V FF {u a +u}' a )^{Lj' a ) 



(61) 



The functions Uff(w) and Vff(^) a re given by the in 
verse Fourier transform of the expressions in Eqs. (|59l 
and can be found numerically. 



B. Corrugation in the second-subharmonic held 
only 

We consider only the forward-propagating fundamen- 
tal held and quasi-phase-matching of the nonlinear in- 
teraction given by the condition § n ^ q + 5^/2 = 0. The 
equations f|55[) can then be written in the following ma- 
trix form 



d 

dz 



A sf (z,t) 
A\ b {z,t) 



= iM 



A Sf {z,t) 
A Sb {z,t) 
At P (z,r) 
AUz,t) 



M(t) = 



r ^72 


K°. 


-UK F 





-K°* 


-S° s /2 








~AiK* F 





-572 


_ K o* 










S° s /2 



(62) 



introducing the nonlinear coupling constant K f (t) 

-A 1 



Kl,a A PF(0,T). 



Eigenvalues Xj of matrix M in Eq. (|d^|) can be derived 
as follows: 



Ai, 2 (t) - V(A0)2 _ 8]Kf1 2 ± 4\KfW4\Kf\ 2 + |^°| 2 , 
A 3 , 4 (r) = -A li2 (T). (63) 

These eigenvalues are real for jif^l < A°/(4y / 2) and de- 
scribe an oscillatory solution. Two real and two pure 
imaginary eigenvalues are found for \K F \ > A°/(4 v / 2) 
reflecting the presence of amplified and attenuated com- 
ponents of the fields. Using the eigenvalues Xj, the solu- 
tion to Eqs. (162[) can be written in a general form using 
operator constants ctj and /3j, 



A Sf (z,t) = ^dj(r)exp[iAj(r)z], 

3=1 
4 

A Sb (z,t) = YsPji^expliXj^z}. (64) 



3=1 



Substitution of the general solution (|M| into the sec- 
ond (or the fourth) equation in (I62p provides the relations 
giving the coefficients j3j in terms of the coefficients &j, 



Pj = - 



K 



0* 



a,, 



J 



(65) 



Assuming real eigenvalues Xj the first (or the third) equa- 
tion in (f62"T) is fulfilled provided that 



«3,4 

7i,2 



71,2^1,21 

AiK F (X 3A + S° s /2) 



\ 2 

A 3,4 



(66) 



The operator constants &i, &2, dj and d| are finally 
determined from the boundary conditions that give us 
the following formulas: 



di(r) 
da(r) 
dl(r) 



MrHr) 



Mi(r) = 
-^/(Ai + ^/^exp^L) 



^ SF (0,r) 
^,(0,r) 
1 1 71 72 

1?1 $2 1?37l ^472 

7i 72 1 1 

^371* nii n 



(67) 
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Using Eqs. (|64l) and (|65)l the solution for the output 
operators can be written as 

4 

A sf (L,t) = ^2aj(T)exp[iXj{T)L], 

4 



l SB (0,r) = 



A" 



-A J (r) + ^/2 



^(r). (68) 



The inverse Fourier transform of the formulas in 
Eq. (1681) and return to the original operators a SF and 
a SB leaves us with the input-output relations written in 
the form: 



E 



-F,B 



du} s U bc (uj s - u s )a 1 s n c {uj s 



+ I dw,V6c(w a +W,)a£t(£^) 



, b = F,B. (69) 



The matrices Ub c and VJ, C in Eq. (|69p can be determined 
numerically in general. 



V. GENERAL NUMERICAL SOLUTION AND 
THE BLOCH-MESSIAH REDUCTION 

To investigate the model numerically we have to re- 
place Eqs. (|8]) for the SSF field by their discrete variants. 
That is why we introduce discrete monochromatic-mode 
operator amplitudes A Sbi that scan the spectral profiles 
with period Alo, 

A Sb>i (z) = y/Au A ab (z,u° s +iAw), 

b = F,B, i = 0,±1±2,.... (70) 

The discrete mode operator amplitudes A ab ^ obey the 
usual boson commutation relations instead of those writ- 
ten in Eq. (j6]). These operator amplitudes can be ordered 
into vectors A SF and A SB . Using these vectors the dis- 
crete form of Eqs. (T5J) can be written as follows: 



d_ 

dz 



A SF (^) 
A SB (z) 



= iM[z) 



M(z) 



A SF (z) 
A\ F (z) 

As B {z) 

M B (z) 

V s -4iJC sF (z) 
MJC* sF (z) -V s 
-JCt 









-ic: 



-V s 4HC sB {z) 







(71) 



Matrix elements of the matrices T> s , K, s , K. s f(z), and 
K, s b{z) are defined as 



V 



s,jk 



Ojk- 



K 



s,jk 



5 jk K s (uj s + kAu 



ICsFjk(z) = AwK n i q (uj° + jAuj, cj" + kAus) 

x ex.p(iS n i tq z)A PF [z, w° + (j + k)Auj], 
IC S Bjk(z) = AujK n i q (uj° + jAuj,uj° + kAu) 



x exp( 



z)A p 



[z,u% + (J + k)Au] 



(72) 



where 5jk means the Kronecker symbol. 

As the matrix A4 in Eq. (|7ip depends on z, only nu- 
merical solution is possible in general. We also need to 
keep quantum features in the solution and so we have to 
solve the system of linear equations (fTTj) for initial vec- 
tors that form a basis. In this way, we reveal the whole 
evolution matrix U that maps the operator fields at z = 
to those at z = L: 



~A SF (L)- 




"A SF (0)" 


Aj F (L) 


= U 


Aj,(0) 


A SB (A) 


A SB (0) 


L M B ( L ) - 




L At s (0) J 



u = 



Uff Ufb 
Ubf Ubb 



(73) 



We note that the operator amplitudes contained in vec- 
tors A SF (L) and A SB (L) obey certain kind of commuta- 
tion relations useful in the numerical solution (for details, 
see [37$ ). 

Partial inversion of the linear relations in Eq. (|73[) re- 
veals the input-output relations among the operator am- 
plitudes, 



"A SF (L)- 




' A SF (0)" 




= U pinv 


A+,(0) 


A SB (0) 


A SB (A) 


L At B (0). 




LAt B (i)J 



U pmv 



Uff —UfbU bb Ubf UfbU bb 
-U bb U B f U bb 



.(74) 



The output operator amplitudes in the vectors A SF (V) 
and A SB (0) obey the boson commutation relations pro- 
vided that the input operator amplitudes given in the 
vectors A SF (0) and A SB (L) fulfil boson commutation re- 
lations. That is why it is convenient to rewrite the re- 
lations in Eq. (1741 into the form of Bogoljubov transfor- 
mation. Using the operators a™ i and a°"* defined in 
analogy with their continuous counterparts we have: 



lout 



u 



V 



linf 



(75) 



The vectors a" 



operator amplitudes a!f, i and a°"*, respectively. 



and a° ut (b = F, B) are composed of the 
and a° ut 

As the matrices U and V describe the Bogoljubov 
transformation, their Bloch-Messiah reduction [3, IH" 
HI can be found, 



U 
V 



XA^Yt, 
XA y Y T , 



(76) 
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where f means the hermitian conjugation and T stands 
for the matrix transposition. The matrices Au and Ay 
are diagonal and contain real nonnegative eigenvalues of 
the decomposition. The matrix Y (X) in Eq. (|76p con- 
tains the right (left) eigenvectors Y^ (X-i). 

The eigenvectors defined by the Bloch-Messiah reduc- 
tion give typical modes of the nonlinear interaction and 
represent a discrete form of eigenmode spectral functions 
(j> n found in Eq. (|35[) in the analytical perturbation ap- 
proach. 

Coefficients B sn and C s „ of the generalized superpo- 
sition of signal and noise [9j written for an n-th eigen- 
mode of the Bloch-Messiah reduction ([76)) and defined 
in Eqs. (j49|) can be expressed using the eigenvalues of 
the decomposition, 



Bs,n 



+ (A[/ jrm Ay jrm C™„_4 + C.C.) , 

» 2 x-iin , »2 x-iin* 
li -U,nn K - / s,nA ~ 1 *-V,nn K - / s,nA 



+ A 



U,nn-L*-V,n7i 



(2B 



s,nA 



i) 



(77) 



Symbol c.c. stands for the complex conjugated term. The 
coefficients B"\ A and C™ nA related to anti- normal oper- 
ator ordering characterize the incident field of eigenmode 
n and can be written as 



B 



s.nA 



a 



cosh 2 (r„) + n„, noisc 
exp(i6>„) sinh(2r„) 



s.nA 



(78) 



where r n is the squeeze parameter, n the squeeze phase, 
and n ni noise the mean number of noisy photons in eigen- 
mode n. The number N s , n of photons in eigenmode n 
[N = (a^a)] attains a simple form: 



^s,n — |£s,n| ~~\~B sni 



(79) 



where £ Si „ gives the initial coherent amplitude in eigen- 
mode n. The principal squeeze variance A Sj „ of eigen- 
mode n is obtained by the formula analogous to that in 
Eq. (O above. 

The importance of the Bloch-Messiah reduction in 
the investigation of squeezing is emphasized by the fact 
that an eigenmode with the lowest value of the princi- 
pal squeeze variance A SFi „ represents the solution of the 
optimization problem for a suitable spectral mode pro- 
file that gives the best possible amount of squeezing (see 
Appendix A). 

In the experiment, either forward- or backward- 
propagating fields are interesting. Their properties can 
be obtained if we decompose the eigenvectors X„ of 
the Bloch-Messiah reduction into their mutually orthog- 
onal forward- (Xf, n ) and backward- (Xs,n) propagat- 
ing parts. The corresponding principal squeeze variances 
\ Sb , n and mean photon numbers N Sb , n are given by the 
weighted sums of the quantities related to the eigenvec- 
tors of the original Bloch-Messiah reduction, 



X, 



N Sb , n = J2° b > nn ' Ns < n '' b = F,B. (80) 

n' 

Using the scalar product • the coefficients Cb jnn > are given 
as 



|xL-x^| 2 

Xft.„ ' X>,n 



b = F,B. (81) 



The number of effectively populated modes belongs to 
the most important characteristics of the SSF field. It 
can be obtained from the analysis of the amplitude cor- 
relation functions (a°^a™ n t ')vac giving the correlations 
between the amplitudes a° u ^ associated with the eigen- 
modes of the Bloch-Messiah reduction. For simplicity, 
the correlation functions are defined for the incident vac- 
uum state | ) vac in the SSF field. Suitability of these 
correlation functions for the determination of the num- 
ber of effective modes originates in the fact that it de- 
scribes paired photons in the SSF field. These paired 
photons are, according to the formula for momentum op- 
erator G in Eq. ([3]), the elementary entities characterizing 
the process of second-subharmonic generation. Using the 
Bloch-Messiah decomposition in Eq. (|76p the amplitude 
correlation functions can be written as: 



( ^nn / X n A\/_ nn Af/_ nn X n . 



(82) 



According to Eq. (|52")l the real number Ay inn A[/ inll de- 
termines the weight of the contribution of an n-th eigen- 
mode. After proper renormalization of these weights 
guaranteeing J2 n (^v,nnA.u,nn) 2 = 1, the number K of 
effectively populated modes is given by the cooperativity 
parameter, 



K 



\Yl n {^V,nnAu,nn) 21 \ 



(83) 



Monochromatic frequency modes play dominant role 
in the experimentally determined quantities. They can 
be easily evaluated using the matrices U and V occurring 
in the solution in Eq. (|T5[) . For example and assuming 
the incident vacuum state in the SSF field, the amplitude 
frequency correlation function N SbSb , tU (ui s j,u) s j') is de- 
termined along the expression 

N SbSh „ u {u} sd ,uj s , f ) = (a° b utt (w SJ )a°" t (w SJ 0)vac 



' d,j'k i 



b,b' = F,B; (84) 



d=F,B k 



uj s .j = ujg+jAuj. The matrices Vbd introduced in Eq. ([84)1 
[and similarly the matrices Ubd used later] are obtained 
from the matrix V [U] in Eq. (|75j) by grouping its matrix 
elements with respect to the propagation direction. The 
spectral photon-number density in mode b is given by the 
quantity N$ bSb>u (u} B j, uj s j ). ^ 

In the time domain, amplitude correlations at two in- 
stants r and r' are characterized by the temporal ampli- 
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tude correlation function N SbS , , t (t, t') expressed as 



N Sb8b ,, T (T,T') = J duj s J duj' s exp(iuj s T)exp(-iuj' s T') 

x a^^KX) 

= (Aw) 2 exp(iw s .jr) exp(— iu> 8 jiT r ) 



(85) 



a) 




For t — t' , the quantity in Eq. ([83)) gives the flux ex- 
pressed in photon numbers. 



VI. DISCUSSION OF THE PULSED 
SQUEEZED-LIGHT GENERATION 

We assume that the incident forward-propagating fun- 
damental field is given by a Gaussian ultrashort pulse 
with the central wavelength A° = 532 x 10~ 9 m and pulse 
duration t p — 1 x 1CP 13 s [see Eq. ([25[l ] originating in 
the second-harmonic frequency generation from a pulsed 
Nd:YAG laser. Knowing its incident power P PF and rep- 
etition rate / (/ = 1 x 10 s s" 1 ), the incident amplitude 
£ p defined in Eq. ([25]) is obtained by the formula 



lP PF Lf3° 



(86) 



b) 




FIG. 5: Intensity spectral profiles |X SjIl | 2 for the scattered (a) 
fundamental and (b) SSF field for the first (plane solid curve), 
second (solid curve with *), third (solid curve with o), and 
fourth (solid curve with A) eigenmode obtained in the model 
with non-dispersion propagation. The profiles are normalized 
such that f duj\X s>n (uj)\ 2 /u)° = 1. In (a) A; = 1.151 x 1CT 7 m, 
A n! = 3.5510 x 10J 6 m and in (b) Aj = 2.459 x 10 _7 _ m, 
A n 



3.5547 x W" m; P v 



1 x l(T a W, L = 1 x 10" 



A. Spectral modes and their properties 



Waveguide's depth t equals 5 x 10~ 7 m. Its width 
Ay is 1 x 10~ 6 m. Depth t\ of the periodic corrugation is 
5 x 10~ 8 m which guarantees single-mode operation at the 
studied frequencies (for more details, see [28[ ). Period A; 
of this corrugation is determined by Eq. ([TUl) such that 
the scattered held a fulfills the resonant condition A° = 
mir/L for the first transmission peak (m = 1), i.e. 6® = 
±^Jn 2 /L 2 + \K®\ 2 . The natural quasi-phase mismatch 
^nTif * s then given in Eq. ([44]) for the fundamental field 
and Eq. (|47[) for the SSF field and determines the period 
A n i of nonlinear periodic poling (for details, see the end 
of Subsec. IIIB). We consider two different waveguide's 
lengths, L = 1 x 10~ 3 m and L = 1 x 10~ 2 m, in the 
discussion. Spectral and modal properties of the SSF 
field are conveniently discussed in the shorter waveguide, 
that provides wider SSF spectra. On the other hand, 
the nonlinear interaction is sufficiently developed in the 
longer waveguide which results in useful values of the 
principal squeeze variance A. 

The section is divided into two parts. In the first part 
spectral modes and their structure are discussed using 
the simplified model with non-dispersion propagation de- 
veloped in Sec. IV. General discussion of the behavior of 
the SSF field as it arises from numerical solution of the 
model of Sec. V is contained in the second part. 



The simplified model with non-dispersion propaga- 
tion is useful namely in revealing spectral properties 
of eigenmodes of the nonlinear interaction. When the 
fundamental-field scattering is considered, the SSF-field 
eigenmodes maintain qualitatively the features obtained 
in the perturbation Gaussian approach in Sec. Ill result- 
ing in the formula ([331) (see also [13, [50]). Thus, an n-th 
eigenmode (n = 1,2,...) has n — 1 zeroes in its intensity 
profile [see Fig. GJa)]. Also, the larger the number n of 
eigenmode, the wider the mode is. 

The eigenmode structure is more complex for the scat- 
tered SSF field that has the forward- and backward- 
propagating components. Intensity profiles of these com- 
ponents in one eigenmode are the same. Moreover, there 
exist pairs of eigenmodes with the same intensity profile. 
However, they are mutually orthogonal due to their dif- 
ferent spectral phase profiles. Also here intensity profiles 
of the components have n — 1 zeros for an n-th pair of 
eigenmodes [see Fig. [5jb)]. 

The absence of inter-mode dispersion leads to the fact 
that spectral widths of eigenmodes are given by the band- 
width of a frequency filter used in the experiment. The 
fact that the propagating monochromatic waves of the 
SSF field are in phase effectively increases the nonlinear 
interaction. As a consequence, smaller values of princi- 
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FIG. 6: Principal squeeze variances A Sjn of the maximum- 
squeezed 15 eigenmodes for the scattered fundamental (solid 
curve with *) and SSF (solid curve with A) field deter- 
mined assuming non-dispersion propagation. Also the princi- 
pal squeeze variances \ SF ,n of the forward-propagating SSF 
field are shown (solid curve with o). For comparison, principal 
squeeze variances A SjIl of the real waveguide without the scat- 
tered fields (solid curve with o) are drawn; P PF = 1 x 10" 6 W, 
L = 1 x 10" 3 m. 



pal squeeze variances A Si „ and greater SSF-field photon 
numbers N StTl compared to the real ones are predicted in 
the model. The principal squeeze variances A s , n of the 
first 15 eigenmodes are drawn in Fig. [5] for the scattered 
fundamental and SSF fields. The principal squeeze vari- 
ances A Si „ for the fundamental-field scattering are smaller 
than those for the SSF-field scattering because the scat- 
tering in the SSF field is weaker [compare the curves in 
Figs. [UJa) and[3Ja)]. In case of the scattered SSF field, 
even the principal squeeze variances X SF _ n characteriz- 
ing the forward-propagating field and given by the for- 
mula (|80p are shown. The values of variances A SF! „ can 
clearly be grouped into pairs which originates in pairing 
of eigenmodes discussed above. However, we note that 
the modes of the forward-propagating SSF field arising 
from the decomposition of eigenmodes into their forward- 
and backward-propagating components are not mutually 
orthogonal. In fact, the number of such modes is twice 
to that given by the dimension of the appropriate space. 
The comparison of values of principal squeeze variances 
A Si „ obtained for the waveguide with and without scatter- 
ing reveals substantial improvement caused by the scat- 
tering (see Fig. [5]). 



B. Pulsed squeezed-light generation 

We analyze the general solution using the model of 
Sec. V assuming the fundamental-field scattering and 
compare the obtained results with those appropriate 
for the waveguide without scattering. Considering the 
shorter waveguide, the intensity spectral profiles of the 
first four eigenmodes are plotted in Fig. [7] for both cases. 
Compared to the profiles of the model without disper- 
sion shown in Fig. [SJa) the obtained intensity spectral 
profiles are naturally bounded by spectral properties of 
the waveguide (material and waveguiding dispersion, see 




FIG. 7: Intensity spectral profiles |X s , n | 2 for the first (solid 
curve without symbols), second (solid curve with *), third 
(solid curve with o), and fourth (solid curve with A) SSF- 
field eigenmode for (a) no scattered field and (b) scattered 
fundamental field. In (a) A„; = 3.5516 x 10~ 6 m and in 
(b) A ; = 1.151 x 1(T 7 m, A ni = 3.5510 x 1(T 6 m; P PF = 
1 x 1(T 6 W, L = 1 X 1(T 3 m. 



Fig. [7]). Whereas the intensity spectral profiles maintain 
the appropriate number of zeros in the case without scat- 
tering, scattering of the fundamental field leads to the re- 
placement of zeros by nonzero minima in these profiles. 
Scattering in the fundamental field considerably broad- 
ens eigenmode spectral profiles [compare Figs. [7f a) and 
(b)] on one side, on the other side it makes the overall 
spectra narrower [see Fig. [9][a) below] . This is caused by 
a complex phase structure of the generated SSF field that 
requires a greater number of eigenmodes in its decompo- 
sition. As these eigenmodes have to be mutually orthog- 
onal, their spectra have to be wider. The mode structure 
does not significantly change when the incident power of 
the fundamental field increases. This is caused by the ap- 
proximation assuming a non-depleted fundamental field. 

On the other hand, increasing values of the inci- 
dent fundamental-field power P PF decrease the principal 
squeeze variances A S! „ in all eigenmodes. This is doc- 
umented in Fig. [5] showing the variances A Sj „ for three 
different values of the power P PF . Comparison of the 
curves in Figs. Uta) and (b) allows to judge the effec- 
tiveness of scattering in the fundamental field from the 
point of view of squeezed-light generation. Scattering of 
the fundamental field not only increases squeezing in the 
eigenmodes, it also considerably increases the number K 
of effectively populated eigenmodes. Whereas K ss 11 
for the case without scattering, K sa 49 is found for the 
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a) 



b) 




n 

FIG. 8: Principal squeeze variances A 3 ,„ depending on the 
mode number n for P PF = 1 X 10~ 7 W (solid curve with o), 
P PF = lx 10" 6 W (solid curve with A), and P PF = lx 10" 5 W 
(solid curve with *) for (a) no scattered field and (b) scattered 
fundamental field. Values of the parameters are the same as 
in the caption to Fig. [JJ 



scattered fundamental field. Larger values of the num- 
ber K of effectively populated modes are important for 
pulsed homodyne detection [7J as they lower the require- 
ments to the amplitude profile of the used local-oscillator 
field. 

As the nonlinear interaction populates a larger num- 
ber of eigenmodes, the intensity spectrum N^ s u of the 
overall SSF field is relatively wide [see Fig. Efa)]. This 
is caused by nearly linear spectral dependencies of the 
linear phase mismatches [see the curves in Figs.[5Jb) and 
in^b)] . On the other hand, amplitude spectral correlations 
given mainly by the fundamental-field spectral width are 
narrow for both considered cases [see Fig. [HJb) for a cut 
across the correlation function N SStU (uj s , wQ]. Also spec- 
tral oscillations originating in dispersion evolution along 
the z axis can be found in these correlations. 

The SSF field is generated in the form of an ultrashort 
pulse (see Fig. [TU]). Whereas its temporal profile is close 
to a rectangular shape for the case without scattering, 
its profile is broken into two parts when scattering of the 
fundamental field is considered. This scattering makes 
the SSF-field intensity spectra narrower and, as conse- 
quence, it also extends the SSF-field duration (from ap- 
prox. 1.3 ps to 2 ps). This extension of field duration is 
also caused by complex spectral phase relations imposed 
by scattering of the fundamental field. Splitting of the 
SSF-field pulse reflects spectral anti-symmetry around 
the central frequencies w° and Period A; of linear 
corrugation and period h. n i,i of nonlinear modulation are 
optimum only for the central frequencies. Whereas their 




FIG. 9: (a) Intensity spectrum iVj sw (w) = N ss ,w(u,U)) and 
(b) cut Ng B IO (u}) = \N S s,w(w®, w)| across the amplitude cor- 
relation function of the overall SSF field for the waveguide 
with the fundamental-field scattering (plane solid curve) and 
without scattering (solid curve with *). The spectra are nor- 
malized such that J dujN SSiU j(uj) /uj° = 1. Values of the pa- 
rameters are the same as in the caption to Fig. [Jj 




FIG. 10: Flux Ns SiT (r) = N ss ,t(t,t) of photon numbers in 
the SSF field for the waveguide with the fundamental-field 
scattering (plane solid curve) and without scattering (solid 
curve with *). It holds that J driV SSjr (T) = 1 and values of 
the parameters are the same as in the caption to Fig. [7] 



values lead to insufficient compensation on one side of 
the spectrum, they overcompensate the nonlinear inter- 
action on the other side of the spectrum. We note that 
the nearly rectangular shape in the case without scat- 
tering is given by the chosen values of fundamental-field 
pulse duration and waveguide length and their relation to 
the group velocities at the fundamental and SSF central 
frequencies. 

In order to reach useful (i.e. sufficiently small) val- 
ues of the principal squeeze variances A Si „, either greater 
fundamental-field powers or longer waveguides have to 
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p p f (io" 6 w) 

FIG. 11: (a) Principal squeeze variances A Si „ and (b) num- 
bers Ns, n of generated photons as they depend on inci- 
dent fundamental-field power P PF . The quantities are shown 
for the waveguide with fundamental-field scattering (A; = 
1.151 x 10~ 7 m, A„i = 3.5510 x 10~ 6 m, solid curves) as 
well as without scattering (A n ; = 3.5516 x 10 -6 m, dashed 
curves). Plane curves are for n = 1 whereas the curves with 
* are for n=23 (with scattering) and n=14 (without scatter- 
ing) equal to the number K of effectively populated modes; 
L = 1 x 10~ 2 m. In (b), the logarithmic y axis is used. 



be considered. Achievable values of the variances X s , n 
as well as numbers N s<n of generated photons depend on 
the fundamental-field power and are plotted in Fig. [11] 
for the waveguide 10 mm long. As shown in Fig. ITT1 
small values of the variances X s . n as well as greater val- 
ues of photon numbers N SiTl are practically found in all 
of the first K eigenmodes. The improvement of squeez- 
ing caused by scattering in the fundamental field is dra- 
matic [see Fig. [TTT a)]. Scattering also increases the num- 
bers N s ^ n of generated SSF-field photons by more than 
two orders in magnitude for the analyzed waveguide [see 
Fig.fTTTb)]. As scattering also increases the number K of 
effectively populated modes (from 14 to 23) the overall 
number N s of generated photons is nearly three orders 
in magnitude greater when the fundamental field is scat- 
tered (see Fig. [nj|. 

Thus, a linear periodic corrugation with suitable pa- 
rameters present in a nonlinear waveguide leads to im- 
portant increase of effective nonlinearity that results in 
great improvement of amplitude squeezing of the gener- 
ated light. 




0.1 i 
0.01 | 

0.001 1 ' • 

0.6 0.5 1.0 1.5 3.0 2.5 

P p f (10~ 6 W) 

FIG. 12: Numbers N a of generated photons depending on in- 
cident fundamental-field power P PF for the waveguide with 
fundamental-field scattering (solid curve) and without scat- 
tering (dashed curve); L = 1 X 10 -2 m. The logarithmic y 
axis is used. 



VII. CONCLUSIONS 

A quantum model of pulsed second-subharmonic gen- 
eration in a nonlinear waveguide with a periodic lin- 
ear corrugation has been developed. Assuming a strong 
fundamental field, the model has been solved for lower 
second-subharmonic field intensities using perturbation 
approach generalized to back-scattered fields. More in- 
tense second-subharmonic fields with negligible inter- 
mode dispersion have been treated by the Fourier- 
transform approach that allows to find partly analyti- 
cal solutions. Numerical approach has been applied in 
the general case. Using the Bloch-Messiah reduction 
spectral eigenmodes suitable for squeezed-light genera- 
tion have been revealed. Scattering by the corrugation is 
more efficient in the fundamental field than in the second- 
subharmonic one. Although scattering by the corruga- 
tion makes the second-subharmonic spectra narrower, it 
broadens the spectral eigenmodes. It also leads to a 
larger number of populated eigenmodes. Phase relations 
in the nonlinear interaction imposed by the corrugation 
also cause splitting of the temporal second-subharmonic 
pulse. In a sufficiently long waveguide, the corrugation 
dramatically increases the number of generated photons 
and, hand in hand, suppresses quantum amplitude fluc- 
tuations. A periodic corrugation thus represents a very 
important and efficient tool for tailoring properties of 
the light generated in modern nonlinear photonic waveg- 
uides. 



Appendix A: An optimum mode for the pulsed 
squeezed light 

We look for a suitable linear combination of the out- 
put operator amplitudes o°"' that minimizes the value of 
principal squeeze variance X s . Using the Bloch-Messiah 
reduction of matrices U and V in Eq. (|76p we can express 
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this combination as follows [6l|: 

<<r X'/ x <;"-- ( A1 ) 

hi 

Coefficients tj fulfil the normalization condition 
Ejl*j| 2 = !• Using Eqs. (gHJ) and g9]) the principal 
squeeze variance with the Lagrange term as a 
function of tj and t* can be expressed in the following 
form: 

A^(t,t*) = l + 2t t A y t-2|t T AyA l/ t|-/it t -t ; 

= 1 + 2t f A v t - 2t f AyAj/t - /it 1 • t. (A2) 

Symbol \x in Eq. (|A2[) denotes a Lagrange multiplier re- 
lated to the normalization of vector t composed of coef- 
ficients tj. 

Derivation of the function in Eq. (| A2|) with respect 
to the coefficients t* and U gives the conditions: 

2 [Ay - Ay Ac/] t = fit, 
2t f [A v - Ay Ay] = fitl (A3) 



Assuming /i = 2(A Vii — Av,uA.u,u) the solution to 
Eqs. (|A3|) is ti = 1 and tj = for j ^ i. As we look 
for the minimum value of X s , we choose i such that its 
principal squeeze variance given by 1 + fi is minimum. 
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